Coulomb blockade in metallic grains at large conductance 
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We study Coulomb blockade effects in the thermodynamic quantities of a weakly disordered 
metallic grain coupled to a metallic lead by a tunneling contact with a large conductance qt- We 
consider the case of broken time-reversal symmetry and obtain expressions for both the ensemble 
amplitude of the Coulomb blockade oscillations of the thermodynamic potential and the 



correlator of its mesoscopic fluctuations for a finite mean level spacing S in the grain. We develop a 
novel method which allows for an exact evaluation of the functional integral arising from disorder 
averaging. The results and the method are applicable in the temperature range S <^ Ec- 



PACS numbers: 73.23Hk, 73.40 Gk, 73.21.La 



I. INTRODUCTION 

The study of electron-electron interactions in meso- 
scopic systems has been at the focus of experimental and 
theoretical interest over the past two decades. One of the 
most striking manifestations of electron interactions at 
low temperatures is the phenomenon of Coulomb block- 
ade.EJ It can be observed by measuring, say the charge of a 
metallic grain which is connected by a tunneling contact 
to a metallic lead and is capacitively coupled to a metallic 
gate that is maintained at voltage Vg, as in Fig. |^. 




FIG. 1. Schematic drawing of a disordered metallic grain 
coupled by a tunneling contact to a metallic lead. The charge 
of the grain is controlled by the gate voltage Vg 

For a poorly conducting contact the charge in the grain 
is quantized at low temperatures and exhibits a charac- 
teristic step- like dependence on the gate voltage Vg. This 
leads to oscillatory gate voltage dependence of all phys- 
ical quantities, referred to as the Coulomb blockade os- 
cillations. The amplitude of these oscillations decreases 
with increasing transparency of the contact and depends 
not only on the total conductance of the contact but on 
the individual transparencies of the tunneling channels. 

In the case of a few transmission channels in the con- 
tact and in the limit of vanishing mean level spacing, 
(5-^0, the Coulomb blockade oscillations of the ther- 
modynamic, transport and thermoelectric quantities are 
completely smeared at perfect transmission.cTa For a fi- 



nite mean level spacing in the grain, the Coulomb block- 
ade oscillations remain finite even at perfectptpansmission 
and exhibit strong mesoscopic fluctuations.H'El 

In the case of a contact with many weakly transmitting 
tunneling channels, the Coulomb blockade oscillations at 
large conductance have been studied in the limit of van- 
ishing mean level spacing inJ:h£. grain, 5, and shown to 
be exponentially suppressedEnld '-^ exp(— 57-/2). Here 
Qt is the conductance of the tunneling contact mea- 
sured in units of e^/(27r?i). In the (5^0 limit for a 
non-random contact the mesoscopic fluctuations in the 
Coulomb blockade oscillations vanish. At finite mean 
level spacing the mesoscopic fluctuations in the thermo- 
dynamic quantities for a weakly transmitting tunneling 
contact were studied in Ref. |l2| and for a multi-channel 
contact near perfect transmission in Ref. |^. In the case 
of a random diffusive contact and vanishing mean level 
spacing in the grain the Coulomb blockade oscillations 
were studied in Refs. 

Here we study the Coulomb blockade oscillations in the 
thermodynamic quantities of a disordered metallic grain 
with finite mean level spacing 5 which is connected to a 
metallic lead by a non-random multi-channel tunneling 
contact and is capacitively coupled to a gate, as in Fig. 
The conductance of the tunneling contact is assumed to 
be large, ^ 1- We assume that electron-electron in- 
teractions in the metallic lead may be neglected due to 
screening and treat the Coulomb interaction of electrons 
in the grain in the framework of the constant interaction 
model: 



Hn = Ec 



(1.1) 



Here Ec is the charging energy, N is the operator of 
the number of electrons in the grain, and q is the di- 
mensionless parameter which is proportional to the gate 
voltage Vg and has the meaning the number of electrons 
that minimizes the electrostatic energy of the grain. We 
consider a disordered metallic grain in which electrons 
move in the presence of a random impurity potential and 
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study both ensemble averaged thermodynamics quanti- 
ties of the grain and their mesoscopic fluctuations. We as- 
sume that the Thouless energy of the grain, Et ~ D/L"^, 
where D is the diffusion coefficient and L is the grain 
size, is greater than the charging energy, Et » Ec- In 
addition, we assume that the time reversal symmetry is 
broken by a magnetic field H and that the cooperon gap 
DeH/{hc) is greater than the charging energy Ec- Un- 
der these conditions the ensemble of disorder potenti^ 
is equivalent to the unitary random matrix ensemble.Ej 

We use the replica formalism to tpaat disorder aver- 
aging in the presence pofjnteractions.ta The application 
of the Itzykson-ZubeiOH integral enables us to evalu- 
ate the resulting functional integral over the Q-matrix 
exactly. The method used here may have useful appli- 
cations to the treatment of non-perturbative interaction 
effects in granulated disordered systems.E3 

The main results of the paper are the analytic expres- 
sions for the avera ge os cillatory part of the grand canoni- 
cal potential, Eq. ( |3.22 ) and for the oscillatory part of the 
correlator of the grand canonic al po tentials at different 
values of the gate voltage, Eq. (3.26). 

The paper is organized as follows: In section || we in- 
troduce the formalism by considering a closed disordered 
metallic grain with non-interacting electrons in the zero 



dimensional limit, and obtain the ensemble average and 
the fiuctuations of the thermodynamic potential f2((?, T). 
In section [II we apply this method to compute the aver- 



age and the mesoscopic fluctuations of the amplitude of 
the Coulomb blockade oscillations in the thermodynamic 
potential n{q,T). We discuss our results in section [V. 



II. THERMODYNAMICS OF 
NON-INTERACTING ELECTRONS 



First, to illustrate the formalism we consider the ther- 
modynamic quantities of an isolated disordered metallic 
grain. In this section we neglect the electron-electron in- 
teractions in the grain. For simplicity, in this section we 
consider spinless electrons. 

We are working within the grand canonical ensemble 
with the chemical potential fi. All thermodynamic quan- 
tities can be obtained from the grand canonical potential 
no{fJ,,T) = -T\nZo{n,T). Here the subscript indi- 
cates that the quantities pertain to non-interacting elec- 
trons and ZQ{fi,T) denotes the partition function which 
can be written as a functional integral over the fermionic 
variables ^ as 



N 



Zo{fi,T) = / ]^ ]^dV'a;,«rfV'a;,nexp 



N 



(2.1) 



In this equation the fermion variables are labeled by 
the Hilbert space indices x and y and by the Matsubara 
indices n and m, the diagonal matrix e = 5xy5nm£n is 
expressed through the fermionic Matsubara frequencies 
£„ = (2n -|- l)7rT, Hq is the single particle Hamiltonian, 
and jjL is the chemical potential. 

We consider the case of broken time-revepsal symme- 
try. Therefore, in the zero dimensional limitEj the single 
particle Hamiltonian Hq can be modeled by a random 
N X N Hermitian matrix belonging to the Gaussian Uni- 
tary Ensemble (GUE) with the probability distribution 



P{Ho) 



exp 



N 

TrH^ 



(2.2) 



As is well knowncEl, the average single particle density of 
states (DOS) i^iE) for such a distribution is given by the 
semi-circle law 



, , N E^ 



(2.3) 



Note we normalized the distribution (2.2) in such a way 
that the width of the semi-circle is equal to unity. Thus, 
throughout this paper all energies are measured in units 



of the width of the semi-circle. For simplicity, we assume 
that the chemical potential n lies near the middle of the 
semi-circle, where the mean level spacing in the dot is 
given by 



N 



(2.4) 



Below we will consider the ensemble averaged thermo- 
dynamic properties p£the dot. For this purpose we resort 
to the replica trickEil and find the averaged replicated 
partition function 



(Zo"(Ai))-(cxp 



t 

T 



(2.5) 



where (. . .) denotes the averaging over the ensemble of 
random matrices Hq wit h t he probability distribution 
function defined in Eq. ( ^.2[ ), and a is the number of 
replicas. The function (Zq (/i)) is initially determined for 
the positive integer number of replicas a and then an- 
alytically continued to a — s- 0. The n-th cumulant of 
the thermodynamic potential ((i7o(/i)")) with respect to 
the distr ibution Eq. (2.2) is then found from (Zq (/i)), 
Eq. (2.5) by using the formula 
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(-l)"da" 

We express the replicated partition function Z" (p) As a result we obtain 



/r, f.\ through the functional in tegr al over the replicated 
fermionic fields tp^ as in Eq. (^.l|) and average Zq (/i) with 
respect to the probability distribution function Eq. (p^). 



a N 



j—l x—l n 



N a 1 ^ 



x,y = l 



n j = l 



(2.7) 



Q to decouple the quartic term in Eq. (2/7) via the 
Next, we introduce a Hermitian a2M x a2M matrix Hubbard- Stratonovich transformation: 



exp 



1 ^ 



2A^ 



x,2/=l 



n j=l 



r In n a 1 

Cax2M / d[Q] exp <^ -yTrQ2 + E E V^x,nQ^m^^,n > 

I x—l n,m i,j — 1 I 



(2.8) 



where at this point to make the sums over the Matsubara fre- 

quencies well defined. Eventually the limit 2M — > oo will 
(2 9) taken to obtain the final expressions. The elements 

of the Q-matrix are labeled by four indices: Ql^jn, two of 

J . 1 ; ■ T-i /R~a\ ■ . 1 u iu i u which refer to the replica space i, j, and two others, n, m 

and the trace m Eq. (E.q) IS taken over both the Matsub- ^ , , jr ,- • 



, \ oaAf(a2M-l) 

Cax2M ~ I ^ ) ' 



J ,. .^r r,7i .r • ,1 1 r refer to the Matsubara spacB. After the integration over 

and the replica mdices. Here 2M is the number of , _ . . . , , . /Ittin • 

1,.^ ; , n . • , ]• TIT 1 -J- c -i the fermionic variables m Eq. E.q) we obtain 

Matsubara frequencies m each replica. We keep it nnite ^ vt_J^ 



ara 



(Zo"(^)) = c„x2A/ / exp \ -yTrQ2 + NTrln (iQ + is + /il 



(2.10) 



where Wje-jntroduced the notation 1 = 5ijSmn. Follow- Eq. ( 2.10 ): Q — Q' — e + ifil. As a result we obtain the 
ing Guhic3 we make the following change of variables in following expression 



-y Tr(Q' - £ + 1^1)2 + NTrln ( iQ' 



(2.11) 



In terms of the new variables U and A the integration 

The matrix Q' can be diagonalized using a unitary ma- measure becomes 

trix U ... 

d[Q']^ A^{A)d[A]d[U], (2.13) 

Q' = U^^AU, (2.12) 
- where d[A] = Yi^Ki d[U] denotes the invariant in- 
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tegration measure in th e un itary group. The Jacobian 
of the transformation ( 2.12 ) is expressed through the 
Vandermonde determinant A(A) = Jl'l-^m ~ ^1)7 where 



the prime indicates that the product is taken over non- 
coinciding pairs of indices {i,rn} 7^ {j,n}. In the vari- 



ables A and d[U] Eq. ( |2.1l| ) acquires the form 



(Zo"(m)) = c„x2M / A2(AMA]d[;7]exp 



-yTr(U-iAU -i + i^it f + NTrhi (^ii) 



(2.14) 



first computed by Harish-ChandrallZl and known ip-the 
Th e in tegration over the unitary group d[U] in physics community as the Itzykson-Zuber integral.tS As 
Eq. ( ^.14 ) can be carried out exactly using the integral a result we obtain 



(^o(m)> 



Cax2M 



aM(a2M-l) 



-NTt - In ( iA 



-(A 

2^ 



i^l)' 



(2.15) 



The Vandermonde determinant A(£) — Jj' {en — £m) 
vanishes if the number of replicas a is greater than unity 
because of the presence of identical Matsubara frequen- 
cies in different replicas. Therefore the expression jn 
Eq. (2.15) appears to diverge. This is not so howeveiES, 



since the integral over d[A] also vanishes when eigenval- 
ues of the matrix e become degenerate. To regularize the 
integral in Eq. ( 2.15| ) we introduce a diagonal regulator 



matrix V = 5ijS„mVi with Vi ^ Vj, and calculate the 
expression 

J 



7r^aM(a2M-l) f ^ J A(A) 



d[K] 



^A(e + V) 
I 



jexp|-A^gF(Ai,,x^j)| 



(2.16) 



where the function F{\lj^,Xm,j) is defined by the equa- 
tion 



,) = - HiX^J -f -(Af„ - x^.,f, (2.17a) 
,j=e^~ili + Vj. (2.17b) 



The ratio of the two Vandermonde determinants in Eq. 
(2.16) is rendere d finit e due to the presence of the regu- 
lators in Eq. ( |2.17b| ): 



A(A) 
A(e-hy) 



n 



XI. 



(2.18) 



{i,m}^{j,n} 



We therefore carry out the integration over d[A] in 
Eq. ( ^.16 ) at finite V and take the limit V" — *■ at the 
end of the calculations. 

For N \ the integration over d[K\ in the r.h.s. 
( ^.16 ) can be performed by the saddle point method.Eja 
To find the saddle point value of at » 1 we may 



ignore the ratio of the Vandermonde determinants in th e 
pre-exponential factor in the right hand side of Eq. (2.16). 
This can be done even after we let the number of Mat- 
subara frequencies 2M go to infinity, as we shall do even- 
tually. We show in Appendix ^ that the shift of the sad- 
dle point equation which arises from the pre-exponential 
factor is small as 5/T and may be neglected. We denote 
the saddle point value of ALby A*„. Minimization of 



the right hand side of Eq. (2.17a) leads to the following 
saddle point equation 



This equation has two solutions for each A* 



Ki = A±(a;,„,,) = — - ± 1/1 



(2.19) 



(2.20) 



Therefo re th ere are 2^*^" saddle points of the integrand 
in Eq. ( ^.16| ). 
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It is clear from Eq. ( ^.16 ) that the fluc tuations of 
about its saddle point value, Eq. ( 2.20 ) are of or- 
der l/\fN. Let us assume for a moment that both the 
temperature T and Vi satisfy the condition <C 
T,Vi — Vj <C 1. Then the ratio of the Vander monde de- 
terminants in the pre-exponential factor in Eq. (2.16) can 
be replaced by its value at the saddle point 



A(A) 



A(£ + V) 



sp 



n 

{i,m}^{j,n} 



(2.21) 



Next, wc formally extend the product in this equation to 
include the diagonal terms | t, m| — {j, n} by replacing 
the ratio in the r.h.s. of Eq. ( ^.21| ) by dX{xm,i)/dxm,i to 
obtain 



A(A) 

A(£ + V) 



= exp <^ - 



1 ^ 



sp 



2 ^ 

irn, 71 j 



In 



— X 



'n,j 



n 

i,7n 



-1/2 



(2.22) 



where the sum over the replica and the Matsubara in- 
dices is unrestricted and includes the terms with m | — 
{j,n}. The factor 1/2 in the exponent in Eq. ( 2.22| ) ap- 
pears due to double counti ng of terms in th e sum mation. 



We now substitute Eq. ( 2.22| ) into Eq. ( 2.16 ) and ex- 
pand the function F{\l^,em,Vi) in Eq. ( 2.17a| ) to sec- 



ond order in the fluctuations around the saddle point 
Eq. (2.20). It is easy to show that 



dXl 



/ dX{Xr 



(2.23) 



Therefore the Gaussian integration over the fluctuations 
around the saddle point, Eq. ( 2.2C ) results in the product 
Y[[{2n/N)dX{x„i,i)/de„i]^^^ which partially cancels the 



last product in Eq. (2.22), and we obtain 



{A} 



exp 

1 •^m.i ) I 



ni.i 



In 



' X{Xr. 



XiXr. 



(2.24) 



where ^ denotes the sum over all the saddle points 

{A} 



given by Eq. ( |2.20D . 

Although, to j ustify the saddle point procedure and 
to arrive at Eq. ( |2.24 ) we as sumed Vi — Vj ^ 1/Vn, 
the saddle point result ( ^.24| ) is exact for the unitary 
ensemble and-iiolds the way to — > 0. It was shown 
by ZirnbauerE^ that this special feature of the unitary 
ensemblejs a consequence of the Duistermaat-Heckman 
theoremr3 

Equation (2.24) is free of divergences in the V 
limit. For the replica-symmetric saddle points it is ob- 
vious because in the "dangerous" terms with m = n the 



argument of the logarithm in Eq. {2-24) remain flnite and 
equal to dX{em — iy^/dsm in the Vi ^ Q limit. For the 
saddle points with broken replica symmetry it is neces- 
sary to sum the contributions of all replica permutations 
of a given saddle point in order to obtain a finite result. 

Next, we take the number 2M of the Matsubara fre- 
quencies to infinity. The resulting infinite sums over the 
Matsubara frequencies in Eq. ( ^.24 ) diverge and need to 
be regularized. The method of regularization can be in- 
ferred from the observation that the sums over the Mat- 
subara frequencies should be regarded as traces of oper- 



ators in the Matsubara space. Let us consider the one 
particle Matsubara Green function as an example 



Gx,y{Tl - T2) = -('0x(Tl)V'y(T2)), 



(2.25) 



where the ave rage is taken over the functional integral 
as in Eq. (2T). It is important to remember that the 

Fermion operators '0 and tp appearing in the Hamilto- 
nian which enter the operator expression for the parti- 
tion function Zo{fj,,T) are normal ordered. Therefore, in 
the time-discretized version of the functional integral for 
the partition function, Eq. (2.1) the Hamiltonian Hq is 
coupled to fermionic fields at different moments of time: 
'ip{ti+i)Hoilj{ti). This forces us to understand the traces 
of the Greens function as 



P N 

TrG= / dr lim Vg,,,(-7?), 



(2.26) 



which coincides with the total number of particles in the 
system. In the frequency representation the regularized 
trace of the Green function is given by 



5 



N 



TrG 



T lim 



+00 

E 

l^ — — OQ X— 1 



5]G,,,(e„)e'^"''. (2.27) 



The exponential factor e*^"'' in Eq. ( 2.27 ) is essential for 
the convergence of the sum over the Ma tsub ara frequen- 
cies. Since the Hamiltonian Hq in Eq. (2J), over which 



the averaging is performed, is normal ordered the traces 
of all the ensuing operators should be regularized by the 
introduction of the factors e*'^"'' as in Eq. (2.27). 



Thus, Eq. ( 2.24 ), where the summations over Mat- 
subara frequencies should be understood according to 
the rule Y^nfi^) ^ I]r=-oo /C*^) *3xp(«'7en), represents 
an exact expression for the averaged replicated partition 
function for N ^ 00. 



Below we will evalu ate E q. ( p.24[ ) at T '3> 6. In this 
case the sum in Eq. (2.24) is dominated by the saddle 
point which corresponds to the lowest value of the ex- 



ponent in Eq. ( ^.24 ). The contributions of all the other 
saddle points are exponentially small in T/d. The leading 
term in the sum in Eq. (2.24) corresponds to the saddle 
point 




(2.28) 



Here the function X'~'{z) = + + is understood 
as an analytic function of z in the complex plane with a 
branch cut from — 2i_to 2i as in Fig. |^, in particular the 
square root in Eq. ( ^.2g ) changes sign when the real part 
of z, diz, crosses zero. 

To obtain the exp ressio n for {Z^{ijl,T)) we take the 
limit 1^ ^ in Eq. (2.24) and perform the summation 
over the replica indices. As a result we obtain 



(Zo"(m)) =cxp<^ lim 



-Na i^(A°(e™-i/i),£m-iM)e'"'"'' 

m— — 00 



-m 



(2.29) 



The first term in the exponent in this equation is linear 
in the number of replicas a and corresponds to the aver- 
age grand canonical potential {D,Q{T,fi)) which is evalu- 
ated in section II A. The second term in the exponent of 
Eq. ( 2.2g| ) is quadratic in a and describes the mesos copic 



fluctuations of ^o{T, jj) which are studied in section II B 



A. Average grand canonical potential 

In this section we evaluate the average grand c ano nical 
potential (rio(/i)). It is easy to show using Eqs. ( ^!^ ) and 
(2.29) that it is given by 



= NT lim V F{\\era - tfJ-),e„, - 

7]-^0 ^ ' 



m— — oo 



(2.30) 



2i+i]x 



Re e 



C 



FIG. 2 
Eq. (|!| 

C in Eq, 



2i+i|J, 

The initial sum over the Matsubara frequencies in 
is re written through the integral over the contour 
2.31). Using the analytic properties of function A" 



in Eq. (2.2i) the integration contour C can be deformed to 



the contour C". 



Th e funct ion i^(A^j£,„ — i/i), s„i~ifi) is defined through 
Eqs. ( ^.17a ) and ( 2.28 ). It has a branch cut from —2i+i^ 
to 2i + in the complex plane of the variable Sm and is 
analytic everywhere else. It is convenient to rewrite this 
sum in terms of the contour integral 



c 



iAi)e*^" ( 



tan f — I 
\2TJ 



(2.31) 
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where the contour C is shown in Fig. |[ Taking into ac- 
count the analytic properties of F{X'^{em — ifJ^), £m — "if^) 



we can deform the integration contour to C and inte- 
grate once by parts to obtain 



NT 
27ri 



de 



In 



(2.32) 



Note that the two terms in the square brackets coincide 
with the trace (in the Hilbert space) of the averaged sin- 
gle particle Green function. The latter is equal to the 
retarded Green function on the right side of the branch 
cut and to the advanced Green function on its left side. 
Therefore, upon the change of variab les e —i{E — fj,) we 
can express the integral in Eq. ( 2.32| ) through th e av erage 
single particle density of states defined in Eq. (^.3[) 



2 

-T J dEv{E) In 



-(£;-m)/t 



(2.33) 



This result is precisely what one expects for the average 
grand canonical potential of non- interacting electrons. 



B. Mesoscopic fluctuations of the grand canonical 
potential 

In this section we consider th e m esosco pic fl uctuations 
of f7o(Ai) at T » (5. From Eqs. and ( |2.29| ) we imme- 
diately find that the second cumulant is deter- 
mined by the second term in the exponent of Eq. ( 2.29| ) 



{{nl{^,)))=T^ 



V" ^i(<^^+^^yn In ( A°(£m-»A^)-A"(£n-z//) 



m.n— — oo 



(2.34) 



r 



As in the previous section it is convenient to rewrite the formula in terms of the contour integral 
summation over the Matsubara frequencies in the last 



mum 



(47ri) 



d£id£2 In 



A"(e,„ - in) ~ \°{en - ili) 



Ci C2 



tan^+i)e*(^i+"^)'' (2.35) 



where the contours of integration Ci and C2 are shown 
in Fig. I 

2i+i|a, 

.C'l 

C2 



±— Re e 



Ci C2 



-2i+i 10, 



FIG. 3. The initial sums over the Matsubara frequencies in 



Eq. (2.34) are rewritten as the contour integrals in Eq. (2.35) 



over the contours Ci and C 2 wh ich are later deformed to the 
contours C[ and C2 in Eq. (b.3f) 



In Eq. ( |2.35| ) the function X°{z - ifi) was defined in 
Eq. ( 2.28D and is understood as a function of z which is 
analytic everywhere in the complex plane except for a 
branch cut from —2i + ifi to 2i + i^x as in Fig. ^. We 
therefore can deform the integration contours Ci and C2 
to C]^, Cj. After that we integrate by parts once with 
respect to both ei and £2 and obtain 
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T2 



{2n) 



deide2Gc{isi + /i, ie2 + /i) In 



-lei/T 



In 



-ie2/T 



(2.36) 



C, Co 



where we introduced the function 



d' 



,2 



Ge(*ei+Ai,z£2+M) = T^ln| ^ ^- ^- |. (2.37) 

deide2 \ £i - £2 



The function G'c(ei,e2) is equal to the smooth part (i.e. 
not containing the oscillations on the scale of the mean 
level spacing 5) of the connected two point correlator of 
Green functions at energies isg + /j, and ie2 + H which was 
obtained by Brezin and ZeeBIl 



Using the fact that the square root in Ao(ei — i/i) has a 
different sign on the different sides of the contour we 
can express ((ilo(/i)^)) through the integral over, say the 
right side of the contours Ci and only. Introducing 
the variables Ei = isi + /i and E2 ~ is2 + /i we obtain 



2 2 



dEidE2Pc{Ei,E2)\n 



'1 + p-(Si-Aj)/T^ 



In 



-(E2-ti)/T^ 



(2.38) 



-2 -2 



where pc{Ei,E2) denotes the sniepth part of the con- 
nected density of states correlatorcZl 



Pc{EuE2 



E1E2 



2n^Ei - E2)^^{A-Ei){A~E. 



(2.39) 



This expression has a singularity when Ei^E 2 which 
leads to the divergence of t he int egral in Eq. ( ^.38 ). If we 
replace pc{Ei,E2) in Eq. (2.38) by the exact connected 
density correlation function pc{Ei, E2) which does not 
diverge at Ei — E2 -^0 the integral in Eq. (2.38) be- 
comes finite and gives the exact expression for the sec- 
ond cumulant of the grand canonical potential. Since we 
have retained only one term corresponding to the low- 
est energy saddle point {Aq} in the sum over t he s ad- 
dle points in Eq. ( |2.24 ), the expression in Eq. ( ^.38 ) is 
missing the oscillatory part of the two point level den- 
sity correlator. These oscil lator y terms arise from the 
other saddle points in Eq. ( 2.24 ) and cut off the diver- 
gence of the integral at the scale of the mean level spacing 
Ei-E2f^6 ^1/N. Therefore the integral in Eq. ( |2.3q ) 
should be cut off at Ei — E2 ^ 1/N which leads to the 
variance ((J7g(/x))) N'^lnN. This result is to be ex- 
pected since the main contribution to the thermodynamic 
potential arises from levels deep below the Fermi surface. 
The fluctuation of the number of such levels is propor- 
tional to IniV and their typical energy is of order —N. 



III. COULOMB INTERACTION AND 
TUNNELING 



In this section we consider a disordered metallic grain 
coupled by a tunneling contact to a clean metallic lead 
and capacitively coupled to a metallic gate as in Fig. [|. 
The Thouless energy of the grain is assumed to be greater 
than the charging energy, Et >> Ec- In addition, 
we assume that the dot is placed in a magnetic field, 
such that the cooperon gap exceeds the charging energy 
DeH/{hc) ^ Ec- It was shown by Efetovlij that under 
these conditions the ensemble of disorder potentials in 
the grain can be described by the unitary random ma- 
trix ensemble. At the same time we assume that the 
Zeeman splitting of electron energy levels, HeH/mc is 
smaller than the temperature T and may be neglected. 
We therefore assume that each orbital level is doubly de- 
generate. 

The single particle levels in the grain are broadened 
due to tunneling into the lead. We assume that the tun- 
neling contact is broad, so that each single particle state 
in the dot is coupled to many lead states. In this case the 
probability distribution P{Tx) of level half-widths F^; is 
sharply peaked about the mean value F. We can there- 
fore neglect the fluctuations of F^: and consider them to 
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be equal to the mean value T which can be expressed 
through the dimensionless conductance of the contact gr 
as 



r = 



9tS 



(3.1) 



We assume that the temperature satisfies the condi- 
tions (5 <C T <C Ec and focus on the Coulomb blockade 
oscillations in the thermodynamic quantities of the grain 
at 3> 1 and at a finite mean level spacing 5. The quan- 
tity of particular interest is the differential capacitance of 
the dot C{q) which can be expressed through the grand 
canonical potential ^{q) of the dot at gate voltage q as 



C{q)^ 



d^n{q) 



(3.2) 



where k is a factor depending on the geometry of the dot. 
In addition to the average differen tial capacitance and 
the grand canonical potential, Eq. ( ^.6| ), we study corre- 
lations of thermodynamic quantities at different values of 
the gate voltage, q and q'. For this purpose we calculate 
the product of replicated partition functions at different 
values of the gate voltage {Z°'{q)Z°' (g')) averaged over 
the ensemble of random Hamiltonians with the distribu- 
The correlators {fl" {q)n"^ {q')) can then be 



tion ( |2.2| ) 
obtained from the formula 



J 



(17"(g)f7™(g')) = (-T)"+™ ( — ) ( — ) (Z"(g)Z" (g')) 



da' 



(3.3) 



A. Average thermodynamic potential 



(n, T) , where /i is the chemical potential of electrons 
in the lead, as a functional integral over the fermionic 
fields tp 



In the absence of the Coulomb interaction but at a fi- 
nite tunneling rate T we can write the partition function 



N oo 



oo N 



/J. * J If 

II II d^l^d^l^exp J2 5] ^ + zrsgn(£)-<5„„(Ho,.y-M'^.y))V^,, 

X — 1 n——oo — - — 1 



n,m——oo x.y—1 a 



(3.4) 



Here the index a denotes the spin of the particle. 

In the presence of the Coulomb interaction, Eq. (1.1), 
the replicated partition function Z"(q) may be written 
as a functional inte gral over the fermionic fields similar 
to that in Eq. (3^). However, in this case the exponent 



acquires a quartic term in the fermion fields ip. A con- 
venient way to proceed is to decouple the interaction via 
the Hubbard-Stratonovich transformation by intrpduc- 
ing an auxiliary field V{t). Following Finkelsteinlij we 
introduce such auxiliary fields Vj in each replica. 



exp ■ 



j Ec{N,{T)-qfdT\^ j d[l/]exp 

^=1 J 




iEc 



* / V,{T)iN,iT)-q)dT 



(3.5) 



where Nj{t) = J2x,a- ''Px,ji''')'^x,ji'^)- Next we denote the 
static part of Vj (r) by Vj and write the integral over Vj 
as a fofn of integrals labeled by a set of winding numbers 

{w}M 



oo 7:T{2Wj+l) 
3 -oo {W^} i T,T(2Wj-l) 

For each set of winding numbers we express Vj (r) as 

V,{t) = V,-2^TW,+^j{t), (3.6) 
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where the phases (/jjiT) satisfy the periodicity condition 
4'jiP) = 0j(O) + 27rVFj. Then, making the gau ge tr an sfor - 
mation '0(r) ^p{T) exp{i4>{T)), using Eqs. (|3.5[), (3.4) 
and averaging over the random matrix Hq with the dis- 



tribution function defined in Eq. (2.2) we write the en- 
semble averaged rephcated partition function (Z^^q)) as 
a sum over the winding numbers 



TTT(2Wi+l) „„ „ / - ~\ f («J -27rrW,)2£iT 

dV ■ r - -fTiQ^+NHln{iQ+'i.j)~j -^^ 



{M/}j = 1^2,(2Vy,-l) 



Here the trace is taken over the spin, the rephca and 
the Matsubara space. The Q-matrix was introduced to 
decouple the quartic interaction of fermions arising from 
the averaging over i?o, denotes the normalized inte- 
gration measure, and we introduced the operator 



J 



where e = i^r, and At,t' = — i/sin[7rr(r — r')]. 

To perform the integration over Q in Eq. (3/?) we pro- 
ceed as in section we make a change of variables 
Q ^ Q' — J and parameterize the matrix Q' in terms 
of the unitary matrix JJ and the diagonal matrix A as 
Q' = U~^AU. The ensuing integral over the unitary ma- 
trices U is of the Itzykson-Zuber type and can be evalu- 
ated exactly. The result is^ expressed through the eigen- 
values J^j of the operator J. The operator J is not diago- 
nal in the Matsubara space but may still be diagonalized 
using a unitary transformation. The index m can be 
identified with the Matsubara index only if = 2-kTW . 
We then perform the integration over A by the saddle 
point method which for the unitary ensemble gives an 
exact result for ^ 1. The saddle point equation for 
the eigenvalues Aj„_^ is given by KnA^m,a - -^m) = 1 
and has two solutions A±(J^) defined in Eq. ( 2.20 ). For 
T ^ S the integration over the static components jii,- can 
be carried out in the saddle point approximation£j The 
saddle point equation reads 



2Ec 



(3.8) 



The left hand side of Eq. (3.8) may be interpreted as the 
number of particles in the dot with the chemical poten- 
tial jJL + i{Vj — 2TTTWj). Note that in order to have a 
reasonably small charging energy the imaginary part of 
the expression in the r.h.s. of this equation should be 
at most of order unity. Since the prefactor in the l.h.s. 
NT = ttT/S is much greater than unity according to our 
assumptions, the saddle point must have an equal num- 
ber of A+ and A_ in order to satisfy (|3.8|). Such a saddle 
point with the lowest free energy is given by 



J* 

•^o(-'^^i) = ^ 



Sgn 



1 



(^m)^ 



(3.9) 



All other saddle points satisfying Eq. (|3.8|) are given per- 



mutations between A+ and A_ in ( B.E ) such that the num- 



bers of A4. and of A_ are equal. These saddle points pro- 
vide only small contributions to the partition function 
~ exp{—T/S) and will be neglected. 



We first solve Eq. (3.8) for Vj at (/)j-(t) = 2TrWjT. In 



this case the operator J is diagonal in the Matsubara 
basis and its eigenvalues are given by 



i/i + Vj - 2ttTWj + Esgn {e„, - 2ttTWj 



P 



Substituting this expression and ( p.9[ ) into Eq. ( |3.8| ) we 
find in the leading order in 6/Ec and 6/T 



V, =2nTWj-i-[q-Noip)], 



(3.10) 



where Nq(ii) is the number of particles in the dot at 
the chemical potential /i given by th e l.h.s. of (3.S) at 



Vj — Wj = 0. To arrive at Eq. (3.10) we have used that 



dNg {fi)/d ^ = 2/(5. It is easy to convince oneself that 
Eq. ( |3.10| ) holds for an arbitrary phase configuration 0j 
in the topological class Wj. We therefore can write the 



operator J at the saddle point ( p. 9]) as 

J = e-i^[q- Noi^i)] + re-'^^(^)A,,,,e'^^(^')(5„-. 

(3.11) 

Retaining only the leading saddle point ( |3.9| ) for each 
set of winding numbers {W} and integrating over the 
gaussian fluctuations of Vj about the saddle point solu- 

tions ( p. 10 ) and using the fact that / (0^ — 2TrTWj)'^dT = 



J 4>']dT + 'in'^TWj we can write the replicated partition 


function as a sum over the winding numbers 



/ X a/2 f ^ J Te^ 

iZ"{q)) - [jk] E / dmf^} e , (3.12) 
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where Z"^^ and A{vf} a-re given by the expressions 
Z{Vj = A{;^}[0]e-^.['"'^^"+^*'*>(*^-^/^], (3.13) 



-Ao(J^) 



(3.14) 



Here we introduced the notation ^{wii^j) for the free 
energy of the dot at the saddle point (3.9) 



J^{i^}(0,) = 7VT^F(Ao(Jii),J^), 



(3.15) 



where F{Xo{Ji), J^) was defined in Eq. ( |2.17a| ). The 
Vandermonde determinant in Eq. (3.14) is generaUy fi- 



nite. This is so because the eigenvalues are non- 
degenerate due to the presence of the auxiliary fields (jfi . 

For the winding number Wj the minimum of the free 
energy -^^j ^f} ('^) achieved on the instanton phase con- 
figuration 02 (r) which can expressed using complex 
variable u = exp(i27rrr) as I 



^ ^ L — It 



(3.16) 



p=i 



Here the complex instanton parameters Zp satisfy the in- 
equahty \zp\ < 1 for Wj > 0, and \zp\ > 1 for Wj < 0. 
The f ree energy J-{w} {4'z) on the instanton configuration 



(|3.16 ) is independent of Zp and is given by 



T 2 ' ^' 



(3.17) 



It is clear from Eqs. (3.17), (3.13 ) th at the largest 
term in the sum over {WYm Eq. (|3.12|) has a mono- 
tonic dependence on the gate voltage, q and corresponds 
to {W} = {0} with aU Wj = 0. The leading oscil- 
latory contributions to this sum arise from the terms 
with {W} = {I}, I — ±1, having only one non-zero 
Wj — I — ±1 and can be chosen in a ways by permu- 
tations between the replicas. Retaining only these terms 
we obtain 



-±1 



/d[0]A{,}[0]exp. 


l^j 1 T 




iEc J j 


/d[0]A{o}[0]exp. 






4Bc J j 



(3.18) 



In Eq. ( 3.18 ) we need to integrate over the fluctu ations of 
(j) around the instanton configuration Eq. (3.16). To this 
end we write (j) = 4'z ~^ 4'j where repr esent s the massive 
modes of the dissipative action in Eq. ( 3.18| ) The integra- 
tion over the zero modes of the dissipative action ( 3.18| ) 
should be performed with the measure (Pz/{1 — \z\'') for 
\z\ < 1, and d2z/[|^|2(|z|2 _ i)] for \z\ > 1. AU the other 
fluctuations (f> have a large mass of order qt- There- 
fore the integration over them may be performed in the 
gaussian approximation. Moreover, in carrying out this 
integ ratio n one may replace the determinants A[(/)] in 
Eq. (3.18) by their values A[02] on the instanton con- 
figurations (3.16) due to their weak dependence on (p. In 



this approximation the gaussian integrals over the mas- 



sive modes in different replicas in Eq. ( ^.Ig ) factorize. 
For the replicas with vanishing winding numbers the in- 
tegrals in the numerator cancel with those in the denomi- 
nator. We are therefore left with the expression involving 
only one replica in which the instanton is present. 



{ZHq)){o} 



= 1 + 2a3?e2''*« 



d^zfiz) A{,}[0, 



Nl<i 



(1 



P) A{o}[0] ' 

(3.19) 



where f{z) denotes the ratio of the integrals over the 
massive modes with and without the instanton 



exp 



2tt^T\z\^ 
'Ecil-\z\') 



J d[(j)\ exp 



/3 {4>fd7 
iEc 



/d[0]exp{-M+;^/^|^} 



(3.20) 
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To arrive at Eqs. ( p^ , ( p^ ) we used that ^ 



2tt T\z\ ^j^g j-^^^ terms with Z = ±1 



Eq. (3.18) are complex conjugates of each other. 

In order to evaluate the ratio of the Vandermonde de- 
terminants in Eq . (|3.19D we need to find the eigenvalues 
of the operator ( |3.11| ) on the instanton configurations. 
This is done in Appendix We then show in Appendix 
^ that in the limit a ^ the ratio Ajij [(/)z]/A{oj. [0] is 
equal to unity. 

The logarithmic divergence of the integral over z for 
short instantons, \z\ 1, in Eq. (3.19) is cut off a t 
1 - ~ T/Ec by the first term iiT/fz), Eq. (|^. 
The ratio of the integrals over the massive modes in 
( |3.20 ) is evaluated in Appendix ^. For very long instan- 
tons, \z\ 0, the result is given by ttit^^tW, whereas 
for short instantons, |z| — s- 1 
QtEc {^ r ^^ ui2 



27r(7rT+r) ' 

and r r the result is 

- ^ (1 - |2p) ln(l + T/ttT)). To evaluate the 
integral ( |3.19D with logarithmic accuracy we may interpo- 
late f{z) for the intermediate instanton lengths between 
the two asymptotics as 



g|.£^cexp|- 



Sc(l-I^P) 



(3.21) 



2n^T{l + ^[l~\z\^]) 
Using this and integrating over z in Eq. ( [3.19| ) we ob- 



tain for the oscillatory part of the ensemble averaged ther- 
modynamic potential 



g'^Ec 



In 



Er 



T + T 



exp ("y) cos(27rg). 

(3.22) 



Eq. ( 3.22 ) is the main result of this section. It is ap- 
plicable in the temperature range i5 <C T <C Ec- For 
r <C T it coincides with the one instanton approxima- 
tion of Ref. |10[ The result of Eq. ( |3.22| ) was obtained 
in the single instanton approximation which holds for 
T:^ g^EceM-^)- 



B. Correlation function (f2(g)n(g')) 

In this section we consider the correlation function 
{il{q)^l{q')) of the thermodynamic potentials at differ- 
ent values q and q' of the gate voltage. For this pur- 
pose we calculate the average product of the replicated 
partition functions at two values of the gate voltage, 
(Z"(g)Z"'(g')). Repeating the steps of section fllA\ 



we 



can write {Z°'{q)Z" (g')) as a sum over the sets of wind- 
ing numbers {W} and {W} 



{Z-{q)Z- (g')) 



EE 

{W}{W'} 



d[(j)]e 



^{V}^{V} nnn 

n,m i—1 j — 1 



7* 



(3.23) 



where and 



are given by Eq. ( 3.13| ). The 
primed quantities in Eq. ( 3.23| ) refer to the set of replicas 
pertaining to the partition function at the gate voltage 
q' . If one replaces the product (Vandermonde determi- 



{W} and {W} in Eq. (3.23) has a monotonic depen- 



nant) in Eq. ( |3.23|) by unity this equation reproduces the 
product of the averaged replicated partition functions. 
The deviations of this product from unity describe cor- 
relations of the partition functions in the different sets of 
replicas. 

For a given set of winding numbers {W} the functional 
integral over the phases is dominated by configurations 
of 4> close to the instanton trajectories, Eq. ( 3.16| ). As in 
the previous section, when integrating over the massive 
fluctuations (j) about the instanton configurations we may 
neglect the d ependence of the Vandermonde determinant 
in Eq. ( 3.23 ) on <f) and evaluate it on the instanton con- 
figuration. The dominant contribution to the sum over 



dence on the gate voltages q and q' and arises from the 
term with all Wj = 0. The largest oscillatory contribu- 
tions to this sum arise from the terms with {W} = {/} 
and {W} — {I'}, 1,1' ~ ±1, having only one non-zero 
Wj = ±1 in each set of replicas. Such terms can be 
chosen in aa' ways. It can be shown, see Appendix 
that the Vandermonde determinant in Eq. (3.23) on 
the instanton configurations with I = I' ^ ±1 goes to 
unity in the limit a, a' —>■ 0. Such terms therefore do 
not contribute to the irreducible correlator of the ther- 
modynamic potentials. The contributions of the terms 
with I = —I' = ±1 are complex conjugates of each other. 
Using this property, denoting the Vandermonde determi- 
nant in the a,a' ^ limit by A;^;/, and using Eq. (3.3) we 
can write the leading oscillatory term in the irreducible 
correlator of thermodynamic potentials as 



{{n{q)n{q')))osc = 23?e2-(«-«') 



d^zd^z'fiz)f{z') / A{i.„i}[^,] 
(l-|zP)(l-|z'P) 1, A{o}[0] 



(3.24) 
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Here f{z) is given by Eq. ( |3.2l|) and we have made an 
inversion transformation for the anti-instanton variable 
z' 1/z'. Therefore the anti-instanton coordinate also 
obeys \z\ < 1. 



The ratio of the Vandermonde determinants ■ 



A{o}[0] 

is determined by the product over the pairs n, m such 
that e„ X e„ < 0, see A ppen dix |^. We are unable 
to evaluate the correlator (3.24) for arbitrary values of 
the instanton parameters and the gate voltage difference 
g — g'.Eil However, for large differences of the gate volt- 
age, q — q' ^ r/(5, the correlator ( 3.24 ) may be evaluated 
with logarithmic accuracy. Indeed, at such gate voltage 
differences even the maximal possible instanton correc- 
tion to th e eigenvalue difference in the denominator in 
Eq. (1^) is relatively small, S{J'^ - 4) ~ T < Siq - q') 
(see Eq. (B7)), and may be evaluated perturbatively. In 



this case 



^{0}101 



only weakly depends on the instan- 



ton parameters z and z' . Therefore the integral is again 
dominated by short instantons 1 — |zp <C T/T for which 
the eigenvalues J/j are accurately described by Eq. 
We show in appendix ^ that for 6{q — q') ^ T,r 



A{o}[0] 



• 4r ~ 



(3.25) 



where b ^ iSjq - q')/ 'iiiT 

Substituting Eqs. ( |3.25[ ) and ( |3.2l[ ) into ( |3.23| ) and in- 
tegrating over z and z' we obtain with logarithmic accu- 
racy the long range asymptotics of the oscillatory part of 
the irreducible correlator of the thermodynamic poten- 
tials 



ler 



5{q - q') 



S{q-q') 
T + ttT 



cos[27r((7 — (?')]. 



(3.26) 



This expression represents the main result of this section. 



IV. SUMMARY AND DISCUSSION 

We studied the effects of Coulomb interaction on the 
statistics of the thermodynamic quantities in an ensem- 
ble of weakly disordered metallic grains with broken 
time reversal symmetry. We assumed that the Thou- 
less energy, Et, in the grain exceeds the charging en- 
ergy, Ec , and considered the case when the grain is con- 
nected to a metallic lead by a tunneling contact with 
a large conductance gx ^ I- We found expressions 
for the oscillatory parts of the average thermodynamic 
potential {il.{q))osc, Eq- (3.22) and of the correlation 
function {rt{ q)n{q ' ))osc at 5{ q 



g') > r,r, Eq. ( p. 261 ). 
The results ( ^.22 ) and (3.26), apply in the t empe rature 
range (5 < E < T. For F < T equation ( |3.22D coin- 
cides with the single-instanton contribution of Rcf. |l^. 
In deriving these results we have neglected the multi- 
instanton contributions. This approximation is valid for 
T g'^Ec exp{—gT/2). The correlator of mesoscopic 
fluctuations, Eq. (|3.26|) , decays as {q~_q')^^\n'^{q — q) 
at S\q - q'\ > T,r. Eqs. ( |3.22| ) and ( |3.26D represent the 
main results of the paper. 

The expressions for the oscillatory parts of the aver- 
age differential capacitance {6C{q)) and of its irreducible 
correl ator (dC{q) 6C{q ')) osc may be easily o bta ined from 
Eqs. (|3.22D and (|3.26|) with the aid of Eq. (^. 



To derive our main results we employed a meth-od. 
based on the extension of the Itzykson-Zuber integraOE^I 
to the infinite dimensional unitary group. In this method 
the Q-matrix degrees of freedom (which describe the 
mesoscopic fluctuations) are integrated out exactly. As 
in the case F ^ 0, at T ^ 5 the replicated partition func- 
tion and its correlator can be written over the topologi- 
cal classes of the phase which deco uples the Coulomb 
interaction term, see Eqs. ( 3.12 ) and ( 3.23 ). Each topo- 
logical class is cKajracterized by the winding numbers in 
different replicasES At large conductance, g ^ 1 the 
functional integral over the phase (j) within each topolog- 
ical class can be taken by the saddle point method. The 
saddle points configurations of the phase ars-character- 
ized by the complex instanton parameters zJ13 At finite 
escape rate from the dot F the functional integral over 
the massive fiuctuations (p around the instanton configu- 
rations depends on the instanton length ~ (1 — |zp 
see Eq. ( |3.21 ). This is in contrast to the F = case 
where it is independent of the instanton length. This 
leads to the fact that the logarithmic divergence of the 
integral over the zero modes of the instantons at low tem- 
peratures is cut off at the escape rate F, see Eq. ( p. 22 ). 
The study of mesoscopic fluctuations of Coulomb block- 
ade oscillations in this approach amounts to evaluating 
the Vandermonde determinants (3.14) on the instanton 
configurations of the phase ( 3.16| ). In the case of the 
average partition function they are equal to unity when 
the number of replicas is taken to zero. In the case of 
the correlator of the thermodynamic potentials the Van- 
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dermonde determinant differs from unity and depends 
on the length of the instantons. The dependence of the 
Vandermonde determinant on the instanton parameters 
z and z' results in the presence of the logarithmic fac- 



tor In^ ^T^j^^T (3-26). In the a-model formulation 



the mesoscopic fluctuations are described by diffusons. 
The dependence of the diffusons on the instanton pa- 
rameters is determined by the z-dcpcndent eig env alues 
of the operator J given by the solutions of Eq. (B6) and 
is easily generalizable-tp time-reversal invariant non-zero- 
dimensional systems£3 The method used here may have 
useful applications to the cr-model treatmeiit|0f interac- 
tion effects in granulated disordered metals.E3 
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In this appendix we demonstrate that the pre- 
exponential factor in Eq. ( 2.16| ) gives rise only to a small 
correction to the saddle point equation, Eq. ( 2.19 ) and 
may be neglected. To show this it is convenient to rewrite 
the integrand in Eq. ( 2.16 ) in the following form 



exp -N^^F{\{Xra.i),Xm,i) + 



In 



(Al) 



Taking a derivative of the function in the exponent in 
Eq. (Al) with respect to A(a;„ji) we obtain the saddle 
point equation for A 



-N 



A(x72^i) ^{^rn.i^ 



0. 



(A2) 



The saddle point value for A, Eq. ( 2.28| ) has been obtained 
under the assumption that the second term in Eq. (A2) 
can be neglected. To check the self-consistency of this as- 
sumption we now evaluate the second term in Eq. ( [A^ ) 



at th e sad dle point A = X^{xm,i), where A° is defined by 
Eq. ( 2.28| ). To this end we rewrite the sum in Eq. ( [A^ ) 
in terms of contour integral 



E 



m— — oo 



— lim / dee''^" f tan f — 



AO(a;„,,0-AO(e)^ 



(A3) 



where the contour C is shown in Fig. and A'^(e) = 
e/2 + ^yT+'e^J4:. Taking into account the analytic prop- 
erties of A(e) we can deform the integration contour to 
C . After a straightforward calculation we obtain 



Y -r- r- = \r; / dx ( tanh — + l) T ^ " In T. (A4) 



A"(£C„,i) 



To obtain the last expression in Eq. (A4) we used the 
fact that A°(a;„,i) - 1 < 1. 

A direct comparison sho ws that the second term in 
Eq. (A_2), given by Eq. (A4), is smaller than the first one 
in the ratio ^ In T (to obtain this estimate we used the 
fact that N = 7r/(5). Therefore we may neglect the sec- 



ond term in Eq. ( |A2|) an d obtain the saddle point value 
for A given by Eq. (|2.28|). 
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APPENDIX B: EIGENVALUES OF J ON THE 
INSTANTON CONFIGURATION 

To find the eigenvalues of the oper ator J it is conve- 
nient to express the operator J, Eq. (3.11), in terms of 
the complex variable u — exp(z27rTr). Introducing the 
notation /i^ = |[g — Nq{^)] we write 



In the z — > limit we obtain that the n = eigenvalue 
is given by Jq (z ^ 0) = Eq — i/ig — T whereas all the other 
eigenvalues are unchanged: J„(z — s- 0) = e„ — i^q + r. In 
the opposite limit of very short instantons, 1— |zp <C TJT 
(as well as for arbitrary z at F <C T ) the sum in Eq. (B6) 
is dominated by a single term in which the denominator 
is small. We then obtain for the eigenvalues with positive 



J + i^lq 

2nT 



-udu 



iT f du' 1 



w u - 



u' 1 



(Bl) 



In Eq. (Bl) we used that the matrices A^.^,-' may be writ- 
ten using the complex variable u as follows: A ~ 



The integral in Eq. (Bl) is taken over the unit circle, 
\u'\ = 1, and the operator should be understood as act- 
ing in the space of functions spanned by the fermionic 
Matsubara frequencies 

For |2;| < 1 the matrix elements of the J in the Mat- 
subara basis are given by 



-2r(l-|znz"z™^?(e„)(?(e,„), 



(B2) 



with d{x) being the step function. It is clear that the 
functions with negative Matsub ara frequencies, 

n < 0, are eigenfunctions of the operator (Bl) with eigen- 
values J„ = e„ — r — i/ig. 

For n > the eigenfunctions of J may be written as 



(u) 



m>0 



(B3) 



where g^"'\u) is a non-singular function of u outside the 
unit circle, \u\ > 1. Using this property and acting with 



the operator J, Eq. (Bl) on (B3) we can write the eigen- 
value equation as 



udu 



r — ifJ,q — Jn{z) 

27rT 



ttT(1-u~^z)9 



(B4) 



Next we introduce the Green's function of the operator 
on the left hand side of this equation. 



G(")(u,«')=E 



m>0 



Jn{z) + i^iq - £„i - r' 



(B5) 



and rewrite Eq. (||) as g{u) = / ^G(w, 
where f{u') is the right hand side of Eq. (B4). Then set- 
ting u — l/z* we obtain the eigenvalue equation in the 
form 



Jniz)^en + T-^^iq-2^\z\''\l 



(B7) 



APPENDIX C: EVALUATION OF THE 
VANDERMONDE DETERMINANT ON THE 
INSTANTON CONFIGURATIONS 

In this appendix we calcu late the Vandermonde deter- 
minant A{vF}[0]i Eq. ( 3.14 ) on the instanton configura- 
tion of the fields 0, Eq. (3.1(:). 

We start with the case of the average thermodynamic 
potential, when the instanton with a unit winding num- 
ber is present only in one replica, j = 1. The correction 
SJn^i to the eigenvalues of J in the j — 1 replica due to 



the presence of the instanton, Eq. (B7) is small in com- 
parison with the band width (which is of order unity in 
our notations). Therefore we may write the correction 
to the Vandermonde determinant Eq. ( p. 14 ), due to the 
presence of the instanton as 



Sin A 




5Ao«i)/9J°i 1 



6Jn 



(CI) 



where Ao(J,° ,,) is given by Eq. (3.9). The terms with 
SmjSn > and with £m,en < in the right hand side of 
this equation vanish. Therefore, to find the dependence 
of the Vandermonde determinant on the instanton pa- 
rameters z we may restrict the product over n and m in 
Eq. (3.14) only to the terms with < 0. 

To evaluate t he lo garithm of the Vandermonde deter- 
minant in Eq. (3.14) in the a ^ limit we can use the 
following considerations. For a given pair of Matsubara 
frequencies n and m in the sum over the replica indices 
i and j there are: i) one term where both eigenvalues 
J* and of the operator J depend on the instanton 
parameter z; ii) 2(a — 1) terms where one eigenvalue de- 
pends on z and the other one does not; iii) (a — 1)^ terms 
where both eigenvalues are independent of z and equal 
to = En + rsgne„ and J^. Thus in the a — + limit 
we obtain 



In A 



-4^ In 



[Uz)-.l^{z)]{J^-J^J 
[J^i-J^{z)][Mz)~-jo^] 



(C2) 



2r|z 



2k 



l-\z 



^ Jniz)+i^iq 



(B6) 



where the summation goes over 7i > 0, m < 0, J,'] = 
En + FsguEn are the eigenvalues of J in the absence of 
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instantons, and J„(z) are given by Eq. (B7). By observ- 
ing that for negative Matsubara frequencies Jm{z) = 
it is easy to see that the right hand side of this equation 
vanishes and Ai[</)z] = 1. 

Next we turn to evaluating the Vandermonde determi- 
nant for the correlator o f ther modynamic potentials, i.e. 
the last product in Eq. ( 3.23 ). In the limit a, a' ^ it 
can be done using the considerations employed above. In 
the case when all the winding numbers Wj have the same 
sign it is not difficult to see that the Vandermonde deter- 
minant is equal to unity. In the case when the instantons 
with both positive and negative winding numbers are 
present this is no longer the case. For the case when one 
instanton with a positive winding number, |z| < 1 and 
one instanton with a negative winding number, |(^| > 1 
are present in the limit a, a' — + we obtain 



A{o}[0] 
^{1 



n 



[Jn{z)-JmmJ^^J?n) 
[J^-J^mJn{z)-J?n] 



(C3) 



For F <C T the eigenvalues J„(z) for \z\ < 1 and n > 
are given by Eq. (B7). The expression for Jm(C) for 
Id > 1 and m < is given by Jm(C) — £m — ^ — 
i^'q + 2F|CP™(ICP — !)■ Substituting these eigenvalues 
into Eq. ( |C3| ) and introducing the new variable z' = 1/C, 
we obtain th e foll owing expression for the exponential 
factor in Eq. (|3.23| ) 



A{o}[0] 

A{1 



-(„ + m + l) + i(,_,') 



4r" 



(C4) 



The sum in the exponent can be evaluated in the lim- 
its when the lengths of the instantons (1 — |zp)/T and 
(1 — |z'p)/r are either shorter or longer than the inverse 
chemical p otent ial difference l/5{q — q'). To evaluate the 
correlator (3.26) with logarithmic accuracy we can i nter- 
polate between the two limits and obtain Eq. (3.25). 



APPENDIX D: INTEGRATION OVER THE 
MASSIVE MODES 

In this appendix we evaluate the ratio of the integrals 
over the massive fluctuations of 6 around the instanton 



configuration in Eq. (3.18). The massive modes may be 
integrated out in the Gaussian approximation. 

We start with the limit of very long instantons, \z\ ^0, 
or |2:| — > oo. In this case it is convenient to work in the 
Fourier basis. We introduce the Fourier components (j)k 
of (/)(t) as 



27rrfc, (Dl) 



and expand the eigenvalues J-j-,-^ of the operator J, 
Eq. (3.11), to second order of perturbation theory in 0^, 



5J, 



(2) 

""^ k^o ^'^ ~ ^^-^ + r(sgnen - sgne„_k) 



(D2) 



Expanding the free energy J^|^j(g, 0) given by 



Eqs. (3.15) and (2.17) in 5J, 



t(2) 



(D3) 



Eq. ( p33| ) we used the fact that 
: at the saddle point. Substituting 



To arrive at 

dF{\j„,j)/d\ -- 

Eq. (p^) into Eq. (D3) we obtain for the second order 



correction to the free energy 



-2iVT^- 



F(sgn£„ - sgnEm) 



(D4) 



In Eq. (D4) the summation over n and m goes over 



n ^ m. Next we collect the terms with k ^ 0. 

To this end we retain only the terms with m — n ± fc in 



the sum over m in the right hand side of Eq. (D4). As a 
result we obtain 



27VrV|0i|2(Ao(Jt 



(D5) 



where we introduced the notation 



1 



27rT/c+r(sgne„-sgn£„_k) 



— 27rT/i:+r(sgnen— sgn£n+k) 



(D6) 

The function Xn,k vanishes for |£„| > 27rTfc. Therefore 
the sum over n in Eq. (D5) is restricted to |£„| < 2TiTk 
and can be straightforwardly evaluated leading to the 
following second order correction to the free energy 



2E 



(27rT)2r|fc| 



k>l 



(5(7rT|fc| +r) 



(D7) 
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Using Eq. (D7) we obtain for the ratio of the integrals over the Fourier components (j)k with fc > 1 in Eq. (3.18), 



c 



ttT + T 



n 

fc=2 



Ec 



irT|fe|+r 



Ec 



77TgT|fc|(|fc|-l) 
7rT|fc|+r 



9tEc 
27r(7rT + T) ' 



(D8) 



In the hmit F ^ T, Eq. (DS) reproduces the resuh of 
Rcf. |l^. 

Next we evaluate the ratio of t he in tegrals over the 
massive fluctuations of (p in Eq. ( |3.18 ) in the limit of 
short instantons, 1 — |zp ^ T/F. To this end we use 



Eq. (D3) to rewrite the correction to the free energy in 
the form 



6F = -7Vr^sgnen(5Jn, 



(D9) 



where SJn up to second order in is given by 



SJn 



(kn|(/'|</5ri 



(DIO) 



here \ipn) are the eigenf unctions of the operator J and 
J„(z) are its eigenvalues. F or th e positive frequencies e„ 
the latter are given by Eq. (B7). 



For a finite length of the instanton the operator J 
is no longer diagonal in the Fourier basis. However, 
in the limit of short instantons, 1 — |zp ^ T/T, the 
non-diagonal matrix elements of J, Eq. ( p32| ), are small 
and may be treated perturbatively. Treating the instan- 
ton correction to J given by the last term in Eq. (p32|), 
7,„,„ - -2F(z*)"z"(l - |z|2)0(e„)0(e„) as a perturba- 
tion we can write the eigenfunctions in Eq. (DlC) as 



J2 



(Dll) 



here = £„ -|- FsgnCn and \n) are the eigenvalues and 
the eigenfunctions of the opera tor J in the abs ence o f the 
instanton. Substituting Eq. (Dll) into Eq. (DIO) and 
expanding the denominator Eq. ( DlC ) to linear order in 
7„.„ after a lengthy but straightforward calculation we 
obtain the correction to the free energy in Eq. (D9) in 
the form SF = X);,;' (l)i(l)i'SFi^i> , where 



SFu, = 



n,i' 



-fF/TTT (1^1 -i-r/^r)(|/'| -t-F/^T) 



(D12) 



Since the functional integral over the fluctuations de- 
pends only on the absolute value of z we as sume d z to be 
real without loss of generality. Equation ( D12 ) is writ- 
ten in the Matsubara basis. It is more convenient to use 
the basis of Ref. ^ which diagonalizes the matrix SF, 
Eq. ( D12| ), at F = 0. Denoting the functions in this basis 
by I (^m), where 



m — 1 
TO > 2 
TO < -2 

TO = —1, 



(D13) 



we can express the matrix elements of SF in this basis as 

6Fnm = J2'^ipn\l)SFi^i,{l'\^m), (D14) 



where {l\ipm) are the elements of the unitary transforma- 
tion matrix between the two bases. 




I > TO 

= TO-1 (D15) 

TO = 1. 



Substituting Eqs. ( |D15| ), (^T|) into Eq. ( |D14| ) after 
straightforward calculations in the limit of short instan- 
tons, 1 — |zp <C T/r, for the matrix elements 5F„„ in 
the new basis we obtain 



SF, 



nm 5T 



T 



n — 1 + a 



{z 



n+m 



(ri — 1 + a) (to — 1 + a) 



(D16) 



Ll' 



Here n,m > 2 and a = T/ttT. For F = the matrix 
6F in Eq. ( D16D becomes diagonal and independent of 
the instanton size, in agreement with Ref. Eq. At finite F 
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and in the limit of short instant ons, 1 — \z\'^ <C T/F, the 



off-diagonal matrix elements in ( D16|) are small and may 
be neglected for the purpose of evaluat ing th e ratio of the 
two integrals over the fields (j) in Eq. ( 3.20 ). Proceeding 
this way we obtain the following expression for the ratio 
of the two integrals over the massive modes in Eq. ( 3.20| ) 



for 1 - 



<T/r 



2tt^T 



1- — (l-z^)ln( 1 + — 



(D17) 



Note that in the limit z ^ 1 this expression|-does not 
contain F and coincides with the result of Ref.t2l. 
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